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1. Summary 


This report describes developments in the following areas: surface 

area and pore size distribution in electrolyte matrices, electron 
microscopy of electrolyte matrices, surface tension of KOH solutions, 
water transport in fuel cells, and effectiveness factors for fuel cell 
components . 

The surface area and pore size distribution, as measured by 

nitrogen adsorption, have been determined for several electrolyte 

matrices made from asbestos or Fybex (fibrous potassium titanate) , 

and for asbestos matrices impregnated with polyphenylene oxide with and 

without a coating of zirconia in polyphenylene oxide on one side. The 

2 2 

specific surface areas found were: asbestos — 44 m /g; Fybex' — 16 m / g; 

2 

impregnated and coated asbestos — 4 to 8 m /g. Pore size distributions 
suggested that the Astropower bag (an impregnated and zirconia coated 
sample which appears to be somewhat superior for some purposes) had 
a more uniform distribution of pore sizes than any of the other samples. 
However, the most significant observation in this area is the large 
reduction in surface area which results from either impregnation or 
coating. 

Scanning electron microscopy at a magnification of 2000 or 5000X 
suggested that the asbestos fiber bundles in the Astropower bag were more 
completely broken down than those in the other samples. It was also 
noted that the coating on the Astropower bag contained relatively large 
flat plates whereas the solid particles in the other coatings appeared 
to be much smaller and less platelike. 
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A recent paper which reported extensive surface tension data for 
KOH solutions was discovered just prior to preparation of this report. 
These new data correlate well both with our experimental data and those 
reported earlier. Data are correlated over a wide range of KOH 
concentrations and temperatures by a polynomial which utilizes a 5 x 5 

matrix of coefficients, but a computer is required for most calculations. 

♦ 

A simpler correlation is being sought. 

Traces of surface active impurities sharply reduce the surface 
— tension of pure KOH solutions, and it is not clear at present whether 
surface tension data for pure KOH solutions are applicable for fuel 
cells. 

The multicomponent diffusion equation has been solved for water 
transport in the electrolyte matrix for an H^-O^ fuel cell with the 
assumptions of constant temperature and pressure (in the matrix) and 
no convection, Ionic hydration has been taken into account. While 
hydration results in significantly larger fluxes of water and hydroxyl 
ions for a given current density than is the case when hydration is 
ignored, the concentration gradients in the electrolyte matrix turn 
out to be quite small and even large errors in estimation of the 
primary hydration number of hydroxyl ion have relatively little influence 
on the concentration gradients. 

A cell for measuring effectiveness factors of fuel cell components 
is under construction; it is based on a modification of the stagnant 
microelectrode. The differential equations for the cell have been solved 
and a computer program worked out to evaluate effective dif fusivities from 
measurements of the limiting current. 



2 » Surface Area and Pore Size Distribution In Electrolyte 

Matrices — M. C. Lee 

The surface area and pore size distribution were measured by the 
methods described in the Fifteenth Semi annual report for nine electrolyte 
matrix samples supplied by NASA. The samples studied bore the following 
designations: Fybex -1 and -2, Asbestos -1, -2 and -3, 4D0573 Sheet, 

and Bag, Astropower Bag, 11D2172 Sheet. The surface area and the apparent 
pore volume are shown in Table 2.1 and the pore size distribution curves 
in the pore size ranging up to about 20 ran in radius are given in 
Figures 1-9. 

One can also measure the pore size distribution by other methods, 
notably by Mercury penetration. In this method, the volume of mercury 
intruded into the pores of a porous sample is measured at several pressures. 
From these data one can construct pore size distribution curves. It 
has been our experience that the mercury penetration technique appears 
to be suitable for the relatively rigid structures of electrodes, but 
that it leads to erroneous results for electrolyte matrices owing to 
their compressibility. It remains to be seen whether pore size distribution 
curves can be generated by mercury penetration measurements when the 
sample has been inpregnated by PPO and then coated. This question 
will be explored in the near future. 

The surface area and the apparent pore volume tests show significant 
differences between the three different kinds of sample groups. From 
the pore size distribution curves it can be seen that there are several 
preferred pore sizes and these are summarized in Table 2.2 below. 
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Table 2.1 

Surface Area and Pore Volume 
of Electrolyte Matrix Materials 


Surface Area Apparent Pore Volume 



m 2 /g 

ml/g 

Asbestos-1 

44.3 

0.148 

Asbestos-2 

44.4 

0.155 

Asbestos-3 

44 .0 

0.227 

Fybex-1 

15.7 

0.0357 

Fyb ex- 2 

16.7 

0,0666 

4D0573 Sheet 

5.5 

0.0186 

4D0573 Bag. 

7.1 

0.0233 

Astropower Bag 

3.9 

0.0214 

11D2172 Sheet 

7.9 

0,0225 


Table 2.2 

Preferred Pore Sized In Electrolyte Matrix Samples 


Sample 

Pore 

Radius 

, R^_ , nm 

Asbestos-1 

3.1, 

5.5, 

12.0 

Asbestos-2 

1.5, 

2.3, 

10.0 

Asbestos-3 

1.0, 

2.5, 

6.0 

Fyb ex-1 

1.2, 

1.7, 

3.0, 5.0 

Fyb ex- 2 

1.1, 

2.5, 

5.0 

4D0573 Sheet 

2.2, 

3.2, 

9.0 

4D0573 Bag 

1.8, 

10.0, 

15.0 

Astropower Bag 

1.4, 

2.5 


11D2172 Sheet 

3.5, 

5.5, 

10.0 
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Since there exists more than one population of pores, the computer 
program described in the Fifteenth Semi annual report was used in 
analyzing the pore size distribution data. In Table 2.3 the computer 
analysis for the asbestos and Fybex materials is given. 

Table 2.3 


Computer Analysis of Pore Size Distribution 
of Fuel Cell Electrolyte Matrix Materials 



Mean Radius Mean Radius 

Local Maxima & Std. Dev. Local Maxima & Std. Dev. 

1.2 1.1 ± 0.2 1.0 1.1 ± 0.2 

1.7 1.7 ± 0.1 

3.0 2.8 ± 0.6 2.5. 2.4 ± 0.4 

5.0 5.1 ± 0.7 5.1 ‘ 5.3 ± 2.0 

7.7 8.6 ± 2.0 


Asbes tos-1 Asbestos-2 Asbestos-3 


Local 

Mean 

Radius 

Local 

Mean 

Radius 

Local 

Mean Radius 

Maxima 

& Std. 

Dev. 

Maxima 

& Std. 

Dev. 

Maxima 

& Std. Dev. 

0.8 

0.8 

± 

o 

• 

o 

P i 
00 

00 

o 

4 

0.04 



1.0 

1.0 

+ 

0.07 

1.1 

1.1 

4 

0.1 

1.1 

0.9 ± 0.2 

1.4 

1.8 

4 

0.4 

1.5 

1.5 

4 

0.1 

1.5 

1.6 ± 0.1 




■ 

1.9 

2.0 

4 

0.1 







2.4 

2.5 

± 

0.2 

: 2.5 

2.6 ± 0.8 

3.1 

3.0 

4 

0.2 

3.3 

5.3 

4 

2.9 



4. 7 

5.4 

± 

1.0 





8.0 

8.8 ± 4.7 


From Figures 1, 2, and 3 we can see that the pore size distribution 
curves for the three samples of asbestos are very similar. Indeed, 
the differences in the graphs may be more indicative of the reproducibility 
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of the experimental measurements than of real differences between the 
samples . 

Figures 4 and 5 indicate that Fybex forms sheets with similar 
preferred pore sizes to asbestos but that the surface area and cumulative 
pore volume are less than for asbestos because of the larger fiber 
size of Fybex. 

The Astropower bag, the 4D0573 bag and sheet, and the 11D2172 sheet 
are all impregnated with polyphenylene oxide (PPO) and coated on one 
side with a suspension of zirconia (and perhaps other materials) in 
PPO. This results in a sharp reduction of the surface area and apparent 
pore volume (all of these samples have a basic framework of asbestos). 

The pore size distribution curves shown in Figures 6-9 reflect the 
changes in surface area and pore volume, but they' also reflect some 
differences between the samples themselves. The cumulative pore volume 
of the Astropower sample is the highest of the group and the curve in 
Figure 6 reveals less structure than is the case for the other coated 
matrix samples. While the difference between samples is small, it 
does appear that the curve for the Astropower sample is significantly 
different from the others. 

Since the Astropower sample is believed to represent a superior 
electrolyte matrix, it seems worth noting here that the pore size 
distributions in these samples appear to be different in the following 
respects: 1) the Astropower sample appears to have a larger number 

of pores of a given size than do any of the other samples, 2) the 
pore size distribution in the Astropower sample appears to be more 
uniform than that of the other samples as revealed by the smoothness 
of the cumulative pore volume versus pore radius plot. It will be seen 
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in the next section of this report that the asbestos bundles appear to 
be broken down more completely and that the range of particle sizes on 
the coated side is greater for the Astropower sample than is the case 
for the other samples. These observations appear to be consistent with 
those made here. 


x 



3. 


Electron Microscopy of Electrolyte Matrix Samples — M. C. Lee 


As noted earlier the last four of the nine samples listed in 
Table 2.1 have been inpregnated with polyphenylene oxide (PPO) and 
coated on one side with a slurry 'of zirconia in PPO. Scanning electron 
micrographs of both sides of these samples have been made, and some 
of these are shown in Figures 10 and 11. 

In Figure 10 we see scanning electron micrographs of the uncoated 
side of the Astropower bag, the 4D0573 bag, the 4D0573 sheet, and the 
11D2172 sheet all at 5000X magnification. Examination of these electron 
micrographs indicates that the asbestos bundles in the Astropower bag 
are more completely broken up than in the other matrix samples. In 
particular, the proportion of fibers of small radius is much larger in 
the Astropower sample than in the others. 

Scanning electron micrographs of the coated side of these samples 
are shown in Figure 11. Although the magnification is too low to 
reveal as much detail as is desired, the difference between the Astro- 
power sample and the others is very obvious. Many large flat plates 
are on the surface of the Astropower sample but these plates show some 
evidence of being platelike aggregates of very much smaller crystallites. 
The larger plates appear to be 5-10 micrometers in length and perhaps 
one micrometer in thickness. Numerous smaller plates, one to a few 
micrometers in length* appear . In contrast, the other electrolyte 
matrices appear to contain more uniformly sized particles of a 
micrometer or so in size, little evidence of a platelike character, 
and an apparently lower concentration of powder in the coating. 
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Electron micrographs at higher magnification were attempted; 
unfortunately the resolution was poor and little or no information 
could be derived from them. This approach will be explored further in 
an effort to ascertain whether the plates in the Astropower matrix 
are porous, and to determine more of the size, shape, and character 
of the particles in the coating of the. other matrices. 



Distribution and Movement of Electrolyte in Fuel Cells — H. Lathanh 

4 . 1 The Surface Tension of KOH Solutions 

In studying the distribution and movement of electrolyte inside 
porous structures it is useful to know the magnitude of the surface 
tension of the electrolyte and its dependence on solute concentration 
and temperature. The surface tension of KOH solutions at different 
concentrations has been measured at room temperature, and a search of 
the literature has revealed a recent extensive investigation by 
Feldkamp (1) for a wide range of temperatures and KOH concentrations, 
with which our experimental measurements agree well at the concentrations 
of interest. 

In this section of the report we shall discuss our own measurements, 
and then we shall deal with those of Feldkamp and other investigators. 

4.1.1 Experimental Methods — The Wilhelmy Plate Method 


Measurements were made on KOH solution in laboratory atmosphere 
by the Wilhelmy plate method, using a surface tensiometer. The method 
consists of measuring the pull, W, of a liquid solution on a completely 
wet table vertical plate of perimeter, p, and the surface tension is 


calculated from 


dynes W dynes 

cm p cm 


(4.1-1) 


This method is only moderately accurate owing to the difficulty 
in positioning the plate very precisely so as to eliminate buoyancy 
forces, and to the presence of a supermeniscus film whose thickness and 
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density depend on the solute concentration. The force W is measured 
by a torsion spring, but it was preferable in our case to calibrate 
the apparatus using pure water for which the surface tension is known 
with high accuracy (2) . 

4.1.2 Preparation of Solutions 

• 

Reagent-grade potassium hydroxide pellets were used to prepare 
KOH solutions of concentrations up to saturation. The methods of 
preparation varied. In the most careful runs, the following procedure 
was used: 

Wash five 250 ml glass beakers with soap and water; carefully rinse 
away all the soap. 

Wash the beakers with concentrated nitric acid for a few seconds, 
rinse with de-ionized water, then with distilled water from a quartz still. 

Allow to dry in laboratory atmosphere. 

Repeat the procedure with five petri dishes which will serve as 
covers . 

To the nearest milligram weigh KOH pellets directly into the 250 ml 
beakers to prepare 50 ml of the desired solution concentration. The 
concentrations are chosen to cover a suitable range, including zero 
concentration. 

Add water slowly, with the beaker on the balance, to make the 
desired weight percent concentration within one percent. 

Allow the system to cool in a water bath, until the desired 
temperature is reached with the Petri dish covers on to minimize evaporation. 

Place the beaker on the balance again and add water to compensate 
for the amount evaporated, to the nearest centigram. 
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Cover the beakers, leave them in the bath again, and prepare the 
surface tensiometer by boiling the hanging plate in distilled water 
and then drying over a gas flame. 

4.1.3 Experimental Results 

Measurements at 21 to 23°C on solutions ranging in concentration 
from 0.5 to 51 w/o yielded a total of 46 experimental points, of which 
only 12 were considered good. The surface tension for water was 
readily reproducible and did not differ for distilled or ordinary tap 
water. The surface tension values for KOH solutions, however, showed 
wide fluctuations which characteristically led to the following behavior: 

1) The scatter in the data was larger for higher concentrations. 

2) All values lay below a clearly defined upper boundary. 

3) At the higher concentrations, the readings showed a rapid 
change with time, always decreasing when this occurred. 

These observations point to the presence of some surface active 
agent in the test solutions, whose effect was dramatized by the high 
concentration of electrolyte, considering that high-concentration 
solutions had a propensity to exhibit much depressed surface tension 
values, at times far lower than that of water. When this did not 
happen, the "good" data lay on the upper boundary, and were considered 
to represent the true values of the surface tension of KOH. These 
data are listed in Table 4.1. Although the precision of the measurements 
was about 1%, the accuracy is lower for the reasons mentioned. Figure 12 
shows a comparison of the experimental data with values reported by 
Faust (3) at 18°C and the correlation of Feldkamp (1) at 20°C, all set 


on a scale relative to water. 



Relative Surface Tension 





25 


Table 4.1 

Surface Tension of.KOH Solutions 
T = 21 - 2 3°C 


Concentration 


o/a 

o 

w/o 

a dynes /cm 


0 

72.5 

1.000 

2.2 

73.8 

1.018 

5.0 

74.1 

1.022 

5.3 

74.9 

1.033 

10.0 

77.0 

1.062 

10.2 

77.1 

1.064 

15.0 

79.3 

1.094 

19.8 

82.1 

1.132 

20.0 

82.7 

1.141 

23.0 

85.4 

1.178 

25.2 

85.5 

1.179 

30.3 

90.1 

1.242 

40.0 

99.3 

1,369 

At the lower concentrations, the experimental data are higher than 

Faust's by 1% and the 

curve from Feldkamp 

is lower by 2%. However, 

for concentrations above 25 w/o, the curves show good agreement. Data 

from the International 

Critical Tables (2) 

along with some others 

between 15°C and 18°C 1 

were plotted in Figure 15 of our previous report 

These agree with Faust 

's at concentrations 

below 10 w/o but diverge at 



26 

higher KOH concentrations. Their range is, however, limited to 17 w/o 
making them of limited importance for fuel cell electrolytes. The 
three sets of values shown in Figure 12 have each been correlated with 
a second-degree equation in concentration: 

Olo = EL + B_m + B 0 m 2 
o 1 2 3 

where: o = surface tension of solution, dynes/cm 

Cf Q = surface tension of water at the same temperature, dynes/cm 
m - KOH concentration, w/o 

Values of the coefficients and the closeness of fit are listed in 
Table 4.2. 


Table 4.2-a 

Coefficients of the Correlations: 



0 = B 1 + B^m 

+ B 3 m 2 



Experimental 

Faust 

Feldkamp 

B 1 

1.0035 

0.9995 

0.9988 

B 2 

0.004406 

0.003135 

0.00231 

B 3 

0.0001170 

0.0001408 

0.0001768 

Goodness of 
Fit 

0.40% 

0.08% 

0.53% 

Concentration 

Range 

40 w/o 

28 w/o 

60 w/o 

Temperature 

21-23°C 

18°C 

20°C 
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Table 4.2-b 

Values Given by the Correlations 


m, w/o 

Experimental 

Faust 

Feldkamp 

Agreement 

20 

1.138 

1.119 

1.116 

2.0% 

30 

1.241 

1.220 

1.227 

1.7% 

40 

1.367 

1.350 

1.374 

1.8% 

50 

1.516 

1.508 

1.556 

3.2% 


4.1.4 The Data of Feldkamp 


Feldkamp (1) has made the most extensive measurements of the 
surface tension of KOH solutions for concentrations up to 60 w/o and 
temperatures from 20°C to 170°C, using the method of maximum bubble 
pressure. His data, which we feel are worth reproducing here, were 
correlated to within 1% by the following equation: 


a(t, O 


5 f 5 

l l 

i=l ^=1 


a.. 

xk 


k-1 


.i-1 


(4.1-2) 


where: t = temperature, °C 

£ = concentration, weight ratio of solute to solution 
and a., is the matrix of coefficients given by Table 4.3. 



Table 4.3 


Matrix of Coefficients of the Feldkamp Equation 


hQ| 

1 

2 

3 

4 

5 

i 

75.4787 

-0.138489 

-0.336392 x 10~ 3 

0.475362 x lO - * 6 

-0.264479 x 10" 9 

2 

-32.8890 

1.34382 

-0.910138 x 10” 2 

0.396124 x 10“ 4 

-0.573565 x 10“ 7 

3 

614.527 

12.8736 

0.104855 

-0.449076 x 10“ 3 

0.651193 x 10~ 6 

4 

-1455.06 

39.8511 

-0.344234 

0.144383 x 10~ 2 

-0.207599 x 10 _5 

5 

1333.62 

-38.3316 • 

0.335129 

-0.137313 x 10~ 2 

0.194911 x 10~ 5 


co 

oo 
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Feldkamp also gave an expression correlating the relative surface 
tension, a/c , with the reduced temperature and concentration, but 
the fit is poor at the lower temperatures. 

Equation 4.1-2 is of the form: 

a(t, 5) = A x + a 2 5 + a 3 5 2 + A 4 £ 3 + A^ 4 (4.1-2) 

where the coefficients are in turn fourth-degree polynomials in 
temperature. In spite of the complexity of this expression and its 
impracticality for hand calculations, it is not readily reducible 
owing to the equal importance of all terms at high temperatures. For 
convenience, Table 4.4 gives values of the coefficients of Equation 
4.1-2 for different temperatures. Table 4.5 gives values of the surface 
tension which are plotted in Figures 13 and 14, emphasizing the concentration 
and temperature dependence, respectively. Figures 15 and 16 show in 
the same respect the behavior of the relative surface tension. 



Table 4.4 


Temperature, * 

>C A 1 

Coefficients of 
A 2 

the Feldkamp Equation 
A 3 

\ 

A 5 

10 

74.061 

-20. 322 

495.83 

-1089.6 

982.46 

20 

72.578 

-9.3454 

395.51 

-789.51 

690.37 

30 

71.034 

0.25742 

311.09 

-532.04 

449 . 79 

40 

69.431 

8.6900 

240.28 

-324.70 

253.67 

50 

67.771 

16.142 

180.92 

-155.60 

95.404 

60 

66.058 

22.788 

131.03 

-18.274 

-31.147 

70 

64.293 

28.792 

88.766 

93.159 

-131.65 

80 

62.479 

34.300 

52.457 

184. 14 

-211.29 

90 

60.619 

39.448 

20.577 

259.59 

-274.81 

100 

58.715, 

44.356 

-8.2401 

323.94 

-326.47 

110 

56.769 

49.131 

-35.203 

381.12 

-370.06 

120 

54.783 

53.866 

-61.365 

434.56 

-408.92 

130 

52.759 

58.641 

-87.625 

487.20 

-445.89 

140 

50.700 

63.521 

-114.72 

541.46 

-483.37 

150 

48.607 

68.558 

-143.24 

599.29 

-523.29 

160 

46.483 

73.790 

-173.61 

662.13 

-567.10 

170 

44.328 

79.242 

-206.10 

730.91 

-615.79 


Table 4.5 


Temperature 


Surface Tension of KOH Data From Feldkamp 

Concentration , w/o 



°C 

0 

10 

20 

30 

40 

50 

60 

10 

74.06 

76.00 

82.69 

91.13 

100.7 

113.1 

132.4 

20 

72.58 

74.88 

81.36 

89.78 

99.59 

111,9 

129.4 

30 

71.03 

73.68 

80.00 

88.39 

98.38 

110.5 

126.6 

40 

69.43 

72.40 

78.59 

86.95 

97.06 

109.1 

123.9 

50 

67.77 

71.05 

77.14 

85.47 

95.66 

107.6 

121.3 

60 

66.06 

69.63 

75.66 

83.94 

94.17 

106.0 

118.9 

70 

64.29 

68.14 

74.14 

82.37 

92.60 

104.3 

116.6 

80 

62.48 

66 . 60 

72.57 

80.75 

90.97 

102.6 

114.3 

90 

60.62 

65.01 

70.97 

79.09 

89.27 

100.8 

112.2 

100 

58.71 

63.36 

69.33 

77.38 

87.51 

98.92 

110.0 

110 

56.77 

61.67 

67.64 

75.63 

85.71 

97.04 

107.9 

120 | 

54.78 

59.95 

65.92 

73.84 ■ 

83.85 

95.14 

105.9 

130 

52.76 

58.19 

64.17 

72.08 

81.96 

93.21 

103.8 

140 

50.70 

56.40 

62.37 

70.14 

80.03 

91.25 

101.8 

150 

48.61 

54.58 

60.55 

68,22 

78.07 

89.28 

99.80 

160 

46.48 

52.73 

58.69 

66.28 

76.08 

87.30 

97.78 

1 7U 

| 44.33 

50.86 

56. /9 

64 . 30 

74.06 

83.80 

95.75 


Surface Tension, dynes/cm 





Surface Tension, dynes/cm 


120 



Temperature, °C . 

Figure 14? Surface Tension of KOH — Temperature Dep endence 
(Data From Feldkamp) 



Relative Surface Tension 


34 



Weight-Percent Concentration 
Figure 15: Relative Surface Tension of KOH- 


Concen t ra tion Dependence (Data From Feldkarap) 
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4.1.5 Discussion 

The agreement between the experimental data and those of Faust and 
Feldkamp bears out the reliability of the former. It is important to 
realize, however, that these represent the ideal situation of pure KOH 
solutions meticulously prepared and measured in an inert atmosphere (1). 
Our observations tend to show that, in many instances, the measured 
values of surface tension are much lower than the expected values, and 
that the difference is greater at the higher concentrations. Since 
the test solutions were prepared carefully, it is likely that the surface 
active agents were present in the laboratory atmosphere and- slowly 
adsorbed onto them. If such is the case, then either the presence of 
ions in the solutions greatly enhances the effect of these agents, or 
it facilitates their adsorption, or both. Until the agents and effects 
are identified and quantified, there is little point in knowing precisely 
the surface tension of KOH under ideal conditions. It is likely that 
in mass-produced fuel cells operating for long times, the organic 
structural materials pollute to a certain degree the electrolyte solution, 
and that this alters its surface tension very significantly. It is thus 
important to evaluate the effect of these pollutants and that of carbon 
dioxide on the surface tension. It is also important to determine the 
dependence on surface tension of the overall cell performance. 

4 . 2 Modelling of Wetting and Penetration of Porous, Hete r ogeneous Media 

Literature and analytical studies of wetting in several simple porous 
or heterogeneous media are in progress but these have not yet led to 
definitive results. They will be described in a later report. 



5. Water Transport In Alkaline Fuel Cells — S. P. Lui 

5 . Introduction 

Alkaline fuel cells, as pointed out by Kordesch (5), are technologically 
the most advanced fuel cell type today. Despite their high power output 
and increasingly lower cost, many problems remain to be solved before 
they can have extensive commercial applications. The concentration 
gradients of electrolyte in the matrix of an operating fuel cell have 
been measured by Miller and Fornasar (6), and by Lundquist and Vogel (7). 
However, to the best of the author's knowledge, the transport of free 
water within the electrolyte matrix has not been discussed in the 
literature. The principal objective of the present work is to look 
into the problem of water transport in the electrolyte matrix and to 
develop ways of estimating the concentration gradients of free x^ater 
and the electrolyte in the matrix. 

The overall reaction product in the alkaline ce H 

water. Water generated by this reaction must be removed from the cell 
to avoid dilution of the electrolyte, but x^ater transport x^ithin the 
fuel cell is also very important in another sense; it is both a 
reactant at the cathode and a solvator of the hydroxyl ions formed 
there. The reactions at the electrodes can be written as follows: 

Anode reaction : 2H ? + 40H . s H 2 0 - (4 + 4s_)H 2 0 

Cathode reaction : 0^ + ^ + 4s_)H 2 0 = 40H~. s_H 2 0 

Overall reaction : 211^ + 0^ = 2H^0 
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where s_ is the primary hydration number of the hydroxyl ion. 

It can be seen, then, that water must be transferred to the 
cathode both to maintain the cathode reaction and to solvate the 
hydroxyl ions produced there. It seems important to note here that 
water transport from the anode to the cathode involves transport of 
solvent water whereas water transport from the cathode to the anode 
involves transport of hydrated hydroxyl ions. 

All of the ions in the electrolyte are hydrated but our principal 
concern is with the hydroxyl ions as these carry most of the intra- 
cell current and they participate in the electrode reactions, whereas 
potassium ions do not. As pointed out by Bockris (8), water may be 
regarded as being bound to ions in two ways : (1) strongly bound and 

(2) weakly bound. Strongly bound water— called the primary hydration 
sheath — is held so strongly to the ion that it moves with the ion. 
Weakly bound water — called the secondary hydration sheath— is attached 
so loosely to the ions that it behaves like a diffuse cloud around the 
ion and is in continuous interchange with water in the surrounding 
medium. We shall consider only the primary hydration sheath here. The 
number of water molecules in the primary hydration sheath of an ion is 
called the primary hydration number of the ion. 

5 . 2 Mass Transfer in the Electrolyte Matrix 

As we have seen earlier, transport processes in fuel cells are 
very important. In the follox^ing treatment mass transfer in the 
electrolyte matrix is considered. The electrolyte matrix usually 
consists of a thin sheet of specially processed asbestos with the voids 
between the asbestos fibers being filled with electrolyte, which, in our 
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case, is concentrated potassium hydroxide. We are able to deal with 
transport in this fuel cell component in a fairly rigorous manner 
because no chemical reactions occur in the electrolyte matrix. All 
transport takes place in the electrolyte solution which has three 
constituents: water, potassium, and hydroxyl ions. However, these ions 

do not exist as bare ions; rather, they exist as solvated ions and the 
water in the primary hydration sheath around the ion is fairly tightly 
bound so that it moves with the ion on the average. The equation of 
continuity for each species in the solution is given by 


3c . 

7r — ■ = -V*N . + R. 

dt xi 


(5.2-1) 


2 

where th is the flux of species i, expressed in moles/cm -sec, and FL 

is the rate of production of species i by chemical reaction expressed 
3 

in moles/cm -sec. Since no chemical reactions occur in the electrolyte 
matrix, is zero, and since the electrolyte matrix is normally very 
thin, we need consider one-dimensional transport only. Therefore, at 
steady state Equation 5.2-1 becomes 


or 



= 0 


N. = constant 

l 


(5.2-2) 


In a multicomponent system, the momentum flux depends only upon 
the velocity gradients and the energy flux depends upon the temperature 
gradient. In the electrolyte matrix considered, although there may be 
some local variations in temperature, the temperature differences across 
the matrix are small compared to the absolute temperature. Therefore 
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we shall assume constant temperature in the electrolyte matrix. Also, 
because the electrolyte is contained in the small void spaces of the 
electrolyte matrix, we shall assume convection to be absent, i.e., the 
bulk velocity of the fluid in the matrix is zero. Hence it is not 
necessary to consider the equations of motion and energy. 

The mass or molar flux of a species depends both on mechanical 
driving forces (ordinary, pressure, and forced diffusion) and on 
thermal diffusion effects. The multicomponent diffusion equation is (9) 


c . 
1 


Vu. + s .VT - -- Vp 
xi P 



(5.2-3) 


where p = electrochemical potential of species i, J/mole; 

£ = partial entropy of species i, J/mole-deg; 

2 

~ diffusion coefficient of i in j, cm /sec; 
ij 

T 

= thermal diffusion coefficient of i, gm/cm-sec; 
c^ = total solution concentration, mole/1; 
v, = average velocity of species i, cm/sec. 


We have already assumed that the temperature in the electrolyte 
matrix is constant and the gradient of pressure is usually very small. 
Thus, the thermal diffusion and the pressure diffusion terms are usually 
of only secondary importance and may be neglected. Equation 5.2-3 
then takes the form 



(5.2-4) 
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Equation 5.2-4 is the generalized form of the S tef an-Maxwell equation (10). 

are the diffusion coefficients describing the interactions 
between species i and j . These coefficients are a set of phenomenological 
coefficients and they follow the so-called "Onsager reciprocal relation”, 
that is, 


= & 
ij Ji 


(5.2-5) 


The term c^Vp^ in Equation 5.2-4 is the driving force per unit volume 
acting on species i. For potassium hydroxide solutions we have two 
independent equations 


C.Vu = RT 

i + 


/ C _L C 

f + O 


; — (v - v, ) + 
C T^o+ 0 + 


c c , 

*“ ~r 


vt. Iv -~ V 


c ?)j = RT 


( c c 

(v ° 


v ) + 


c - c + 

C T*^f 


(V + - V_) 


where 


c 

o 



y 


concentration of total water, moles/1; 
concentration of potassium ion, moles/ 1; 
concentration of hydroxyl ion, moles/ 1; 
total solution concentration, moles/1; 
electrochemical potential of the potassium ion; 
electrochemical potential of the hydroxyl ion. 


(5.2-6) 


(5.2-7) 


The chemical potential of the electrolyte may be expressed as 


h e - v + U + + vji 


- vRTln (cf + a+) 


(5.2-8) 
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where V 

"r 


number of potassium ions produced by the dissociation 
of one molecule of potassium hydroxide; 


V = number of hydroxyl ions produced by the dissociation 

of one molecule of potassium hydroxide; 

V = number of moles of ions into which a mole of electrolyte 

dissociates; 

f = mean molar activity coefficient of the electrolyte; 

0 

a + - a constant of proportionality. 


The number of moles of ions into which a mole of electrolyte dissociates, 
V, is given by 

v = V + v (5.2-9) 


The concentration of the electrolyte is defined by 



(5.2-10) 


The flux equations for the potassium ion and the hydroxyl ion may be 
obtained by rearranging Equations 5.2-6 and 7. For steady state conditions, 
these flux equations can be written as follows: (see Appendix 1 for 

details) . 


N + - 


C V = 

+ + 


V Jr'a c 

" vRTc VV e + 


it 


+ 


+ . 

TT + C ' V 

F +o 


(5.2-11) 


N = c v = - 


V &QL c it 

zrtr V]i + + C V 

VRTc H e z F - o 


(5.2-12) 


where is the diffusion coefficient of the electrolyte and is given by 


A (z , ' z ) 

{£. = _S+ o~ + 


Z - z 

4 * Ot — o - 


(5.2-13) 
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The transference numbers (with respect to the solvent velocity) are 


defined as follows 


z Er 
+ o+ 


t , - 1 - t - jrr — 

+ Z , S' , - Z S' 

+ 0 + - o- 


(5.2-14) 


where we note that z has a negative value and the term -z_ becomes 

<y 

positive. The measured diffusion coefficient, D, is related to by 




D = 


1 + 


dlny + 
din m 


(5.2-15) 


where y+ is the mean molal activity coefficient and m is the molality. 

It can be shown that the relation between the gradient of chemical 
potential and the gradient of concentration is given by (see Appendix 2) 


B'c T c 

— Yp = D 
vRTc e 


1 - 


dine 
c 

dine 


Vc 


(5.2-16) 


Substituting Equation 5.2-16 into Equations 5.2-11 and 12, we get 

it . 


and 


N + - - V + D 


N = - V D 


( dine 1 
1 - 


dine j 
dine 


1 - 


dine 


Vc + — v + c v 
z_j_F + o 


it 

Vc + — ~ + c v 
z F - o 


(5.2-17) 


(5.2-18) 


For water, the total flux consists of the flux due to free water, that 
due to water of hydration, and that due to convection. Since potassium 
ions do not enter into the electrode reactions, the net flux of 
potassium ions is zero and the flux of hydration water with potassium 
ions is, therefore, zero. Thus, we have 

dc 


N H 2 0 ° dx 


fw 


+ s N + c v 
- - o o 


(5.2-19) 
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where is the concentration of free water expressed in moles/1 
and s is the hydration number of the hydroxyl ion. 

We have already assumed that convection is absent in the electrolyte 


matrix and that all transport is essentially one-dimensional, so the 
flux equations become 


N + = -V + D 


1 - 


dine 
c 

dine 


dc 

dx 


z + F 


N = -V D 


dine 
o 

dine 



it_ 

7T. 


(5.2-20) 


( 5 . 2 - 21 ) 


V ■ -» 


dc 

dx 


fw + sN 


(5.2-22) 


In order to solve these equations, it is necessary to know the 
transference number and the hydration number of the ions. In the next 
two sections the transference number and the hydration number in 
concentrated potassium hydroxide solution will be discussed. 

In Section 5.5 the concentration gradients for both the electrolyte 
and free water are calculated using both the transference number and 
the hydration number estimated for 30 w/o KOH. 

5.3 Transference Numbers In KOH Solutions 


The transference number of an ion, t^. 


total current carried by that ion, that is 


t 


i 


i. 

i 



is the fraction of the 


(5.3-1) 


where i. is the portion of current carried by the i species, expressed 

2 

in amp/cm . Also, we note that 
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1 t. = 1 


(5.3-2) 


Since the current density of an ion is proportional to the mobility 
of the ion, it can be shown that the transference numbers of the 
potassium and hydroxyl ions are given by 


t, = r 

+ U, + u 

+ 


and 


t = 


u + + u 


(5.3-3) 


(5.3-4) 


where u^ and u_ are the mobilities of the potassium and hydroxyl ions, 

9 

respectively, expressed in cnT/sec-volt • The mobility of an ion in 
a solution is directly proportional to the conductance of the ion; 
hence transference numbers may also be expressed in the form 

A 


+ 


e + “ 1 - " x. + A 

+ 


(5.3-5) 


where A_^ and A are the ionic conductances of potassium and hydroxyl 
ions, respectively. The equivalent conductance of a binary electrolyte, 
A, is given by 


A = 


Hence , 


Similarly, we have 


A^ + A 


t + ~ A 


t = 


A 


(5.3-6) 


(5.3-7) 


(5.3-8) 


5.3.1 Methods for Measuring Transference Numbers 


The chief methods for the determination of transference number are: 
the Hittorf method, the moving boundary method, and the electromotive force 


method. 



46 


The Hittorf Method 

This method is based on a procedure suggested by Hittorf in 1853 
in which the transference number is determined by measuring the changes 
of concentration in a three-compartment cell in which a given current 

is passed through a solution of the electrolyte. This procedure is 

# 

subject to many difficulties in practice, and accurate results are 
difficult to obtain. However, by using an improved apparatus, Maclnnes 
and Dole (11), and Jones and Bradshaw (12) were able to obtain quite 
accurate results. 


The Moving Boundary Method 


In the moving boundary method, the transference number is obtained 
by measuring the volume swept through by a boundary between two solutions 
during the passage of a given quantity of electricity. The transference 
number of the cation is given by 


* 

FdAc 

t — 

+ Q 


(5.3-9) 


where 


d = 

A = 
* 

c = 


Q - 


distance moved by the boundary, cm; 

2 

cross-sectional area of the tube, cm ; 
concentration in equivalent per cm^; 
quantity of electricity, coulomb. 


For a detailed account of the theory and techniques of the moving 
boundary method, the reader is referred to Maclnnes and Longsworth (13). 
Tliis method gives more accurate results than does the Hittorf method, 
and it is simpler experimentally. 
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The emf Method 

This method of measuring transference numbers is based on measure- 
ments of the emf of concentration cells with and without transference. 

A concentration cell with transference consists of two identical electrodes 
immersed in solutions of the same electrolyte having different compositions 
or solutions of different electrolytes having the same concentration. 

The difference in composition of the two solutions at their junction 
leads to diffusion and the development of a liquid junction potential. 

An example of a concentration cell with transference is 


Ag, AgCl | KC1 || KC1 1 AgCl, Ag 
a l a 2 

and the emf of the cell is given by 


± t 


VRT 

V± Z ±F 



( 5 . 2 - 10 ) 


where the subscript on E indicates that transference takes place, 
and and are the mean activities of the electrolyte in the left 
hand and right hand compartments of the cell. 

A concentration cell without transference is one in which both 
electrodes are the same and immersed in the electrolytic solutions 
of different compositions but with no liquid junction between the 
anolyte and the catholyte. Since there is no liquid junction, ions 
cannot pass from one solution into the other. There is, therefore, 
no liquid junction potential in concentration cells without transference. 
An example of a concentration cell without transference is 
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Ag, AgCl | KCl | K(Hg) | KC1 | AgCl, Ag 
a l a 2 

where it is seen that the connection between the anolyte and the 
catholyte is a potassium amalgam. The emf of the cell is given by 


E « 


+ VRT 
“ V ± Z ± F 



(5.3-11) 


where V is the total number of ions, and V, or V is the number of 

T — 

positive or negative ions produced by the dissociation of one molecule 
of electrolyte; z + or z is the charge number of the positive or 
negative ions. 

Since the difference in the emf of concentration cells with and 
without transference arises out of transference across the liquid 
junction, it follows, then, that we can determine the transference 
number by comparing the emfs of the respective cells. Thus, from 
Equations 5.3-10 and 11 we have 



(5.3-12) 


From Equations 5.3-7 and 8, the transference number of an ion 
is expressed as the ratio of the ionic conductance, A_^, and the 
equivalent conductance of the solution, A, where the values of the 
conductances are those at the particular concentration of interest. 

Since the ionic conductances are concentration-dependent, it is not 
surprising to find that the transference number also vary with the 
concentration of the electrolyte. 

Empirical equations for transference number as a function of 
concentration have been proposed by Jones and Dole (14), Longsworth (15), 
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and Owen (16). These empirical equations give satisfactory results 

for aqueous solutions of 1—1 chlorides but the validity of these 

equations for potassium hydroxide solutions has not been tested adequately. 

According to the Onsager equation the equivalent conductance of 
an ion,X_£, is given by 

A. = A° - (l/2a + BA°) Jl (5.3-13) 

11 1 

where A? is the ionic conductance of species i at infinite dilution, 

l 

2 

expressed in mho- cm /e q, ct and 3 are constants. Thus, the transference 
number of the potassium ion is given by 

x A° - (l/2a + BA°)/c 

fc + = A + A = A - (a + BA (5.3-14) 

+ - O O 

For water as solvent, a = 60.20 and 3 - 0.229. The ionic conductances 

2 

for the potassium and hydroxyl ions are 73.52 and 197.6 mho-cm /eq, 
respectively. Using these values of the constants in Equation 5.3-14, 
values of the transference number for the potassium ion can be calculated. 

Although dilute electrolytic solutions have been studied intensively, 
studies in concentrated electrolytic solutions have been less frequent; 
and while the Onsager theory predicts the concentration depencence of 
the transference number, it is applicable only at low concentrations. 
Attempts to extend the Onsager equation to higher concentrations have 
not met with success. 

Recently, Merenkov (17) measured the transference number of potassium 
ion in KOH solutions from one equivalent to ten equivalents per liter. 

He found that the transference number of potassium ion was independent of 
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temperature from 18° to 65°C but that it decreased with increase in 
KOH concentration. The transference number of potassium ion was 
correlated by the empirical equation 

t, = 0.26 - 0.047 {SZ + 1) (5.3-15) 

-r 

On the other hand, Knob el (18) determined the transference 
number of the potassium ion in KOH solutions by the emf method; 
he found that the transference number of potassium ion was 0.2633 and 
that it was independent of concentration over the concentration range 
from 0.03N to 3. ON. The predictions of the Onsager equation and the 
experimental results of Knob el, and Merenkov are summarized in Table 5.1. 

In fuel cell operations, we are dealing with concentrated 
solutions of KOH. The use of the Onsager equation certainly leads to 
erroneous results. The empirical equation proposed by Merenkov gives 
the concentration dependence of the transference number for concentrated 
solutions, yet it cannot be used for dilute solutions. The value at 
infinite dilution as predicted by Merenkov’ s empirical equation is 
0.213, which is far below the theoretical value of 0.271. At a 
concentration of 1 N, however, the agreement between the Onsager and 
Merenkov equations is relatively good (t^ = 0.176 vs t_^ = 0.166). 

It is clear that Knobel’s data deviate from the Onsager equation very 
significantly when the potassium hydroxide concentration is greater 


than about 0.1 N. 
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Table 5.1 


Summary 

of 

Transference Number 

Data for KOH 

Solutions 



Transference number of 

K + from 

Concentration 
KOH Solutions 

of 

00 

Onsager Eq. 

Knobel Eq. 

Merenkov Eq. 

0 


0.271 

- 

0,213 

0.03 


0.261 

0.2633 

0.205 

0.1 


0.251 

0.2633 

0.198 

0.3 


0.232 

0.2633 

0.187 

1.0 


0.176 

0.2633 

0.166 

3.0 


- 

0.2633 

0.132 

5.0 


- 

- 

0.108 

7.0 


- 

- 

0.088 

9.0 


- 

- 

0.072 

10.0 


— 

- 

0.064 


5.4 The Primary Hydration Number In Concentrated KQH Solutions 

The process of solvation plays an important part in most solution 
phenomena. The hydration number of an ion is the number of water 
molecules attached to the ion, and the absolute value of the hydration 
number has an important bearing on the amount of "free" water present 
in electrolytic solutions. For example, in 30% KOH solution, if the 
hydration number of both the potassium and hydroxyl ions is equal to 
one, then the concentration of "free" water is 36,4 moles/1, while if 


the hydration number for the ions is equal to four, there will not be 
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any ’’free" water available. The concentration of "free" water for 
different values of the hydration number at various potassium hydroxide 
concentrations is tabulated in Table 5.2 (assuming that the hydration 
number of the potassium ion is equal to that of the hydroxyl ion). 

Table 5.2 

9 

Concentration of "Tree" Water 
As A Function Of Hydration Number 


Concentration 
of KOH Solution 
Wt.% moles/1 

Concentration 
s=0 1 

of free 
2 

water , 
3 

c_ , mo 
fw* 

4 

les/1 

5 

10 

2.3 

64.5 

59.9 

55.3 

50.7 

46.1 

41.2 

20 

4.6 

57.3 

48,1 

38.9 

29.7 

20.5 

11.3 

30 

6.9 

50.2 

36.4 

22.6 

8.8 

v- 

r - 

AO 

9.2 

43.0 

24.6 

6.2 

T> 


-• 

50 

11.5 

35.8 

12.8 






See Appendix 3 for sample calculations. 

The concentration of "free" water as a function of KOH concentration 
is shown in Figure 17 for several arbitrary values of the hydration 
number. 

Many workers have investigated the concept of ionic hydration, 
and several experimental methods have been used for determining the 
hydration number of ions. Some of these are discussed below; however, 
it should be noted that the different experimental methods give contra- 
dictory results because of a lack of a really precise definition of 
\ 

the hydration number. As pointed out earlier, solvation is arbitrarily 
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divided into two different categories-primary and secondary. 

Of the several experimental methods for determining primary 
hydration numbers which have been used, the more widely used appear 
to be: ionic entropy, ionic mobility, density, compressibility, and 

the nuclear magnetic resonance method. Some of the experimental 
results which have been reported for the hydration number of the 
potassium ion are summarized in Table 5.3. 

Table 5.3 

+ 

Primary Hydration Numbers of K 
By Various Experimental Methods 

Method Hydration 

Number 


Entropy, Ulich (19) 3 

Density, Darmois (20) 1 

Compressibility , Passynski (21) , 6-7 

Bockris (22) 3.5 

nmr , Hindman- (23) 2.1 

Creekmore et al, (24) 4.6 

X-ray diffraction, Brady (25) 4*0 


In general there is a great variance in the experimental values 
of the hydration number because different methods measure different 
types of solvation (primary or secondary, or a combination of both). 

Of all the techniques available, nuclear magnetic resonance (nmr) and 
compressibility appear to be the most promising methods for determining 
the hydration number of ions. 
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The use of the ninr techniques for determining the hydration 
number of ions is still an active field of research, and several 
approaches can be used to determine the hydration number form nmr 
measurements. Swinehart and Taube (26) determined the solvation 
number of certain diamagnetic cations in methanol from the ratio 
of areas under the nmr absorption curves of the free and solvated 
molecules. It has been shown by Jackson, Lemons, and Taube (27) that 
adding a paramagnetic ion to an aqueous solution containing diamagnetic 
ion produces two separate 0^ nuclear magnetic resonances, one for the 
labile water and the other due to water bound in the primary hydration 
sheath of the diamagnetic ion. The hydration number of the cation is 

then determined either by comparing the areas under the two resonances 

17 . . 

or by measuring the 0 nmr shifts in solutions containing a para- 
magnetic ion and a diamagnetic ion and in a solution containing only 
the paramagnetic ion. The hydration number of ions can also be determined 
by nmr techniques from proton relaxation time measurements. 

Although the nuclear magnetic resonance method has become an 
increasingly important technique for the study of hydration of ions, the 
feasibility of the technique for measuring the hydration numbers of 
potassium and hydroxyl ions in concentrated solutions requires further 
investigation. 

The compressibility method, developed by Passynski (21) also 
appears to be a satisfactory method for determining the hydration 
number. By assuming the compressibility of water in the primary 
hydration sheath to be zero, the sum of the hydration number of ions 
in solution can be determined from compressibility measurements on the 
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solvent and solution. Isemuro and Goto (28) used- the method of 
compressibility to estimate the hydration number of alkaline halides. 

As one might expect, they found that the hydration numbers of 1-1 
electrolytes decrease with an increase in the temperature, but the 
concentration dependence of the hydration number was not investigated. 
Recently, Bockris (22) used the compressibility method to determine 
the sum of the hydration numbers of ions in solutions of halides. By 
combining compressibility measurements with measurements of the ionic 
vibration potential, individual hydration numbers for certain cations 
and anions were estimated. According to his results, the hydration 
number of potassium ion is 3.5 and that for hydroxyl ion can be estimated 
to be equal to 4.3. If these hydration numbers remain constant in 
concentrated solutions, as pointed out earlier, there would not be any 
"free" water available in 30% KOH solution. It is obvious, therefore, 
that the hydration number must vary with the concentration of the 
solution, but there is little experimental information concerning 
the concentration dependence of the hydration numbers of the potassium 
and hydroxyl ions at high concentrations. 

From some of the limited information available in Bockris l s 
paper, the hydration number of the potassium and hydroxyl ions 
may be estimated by extrapolating the available data. This procedure 
results, for 30% KOH solution, in an estimated hydration number of 
potassium ion of 1.5 and 2.0 for the hydroxyl ion. These values are, 
however, very tentative owing to the extrapolation required. Since they 
are the best estimates available, however, they are used in the next section 
to estimate concentration gradients in the electrolyte matrix for 


different current densities. 
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5 . 5 Concentration Gradients in the Electrolyte Matrix 


Having estimated the hydration numbers and the transference 
numbers of the ions, the flux equations for the potassium, ion, the 
hydroxyl ion, and water can be solved. The concentration gradients of 
the "free” water and that of the electrolyte in the electrolyte matrix 
can be estimated as follows: 

For the "free” water concentration gradient, we have from 
Equation 5.2-22 



dc 


-D 


fw 


dx 


+ s N 


i 

F 


Hence the "free” water concentration gradient can be expressed in 
terms of the hydration number of ions, the diffusivity, and the current 
density, that is 


dx ~ DF 


(5.5r-l> 


Values of the diffusion coefficient of the electrolyte at various 
concentrations, taken from Bhatia’s results (29), is shown in Table 5.4. 


Table 5.4 


Diffusivity of KQH Solution at 25 
Concentration of KOH 


wt.% 

moles/1 

D x 10 

10 

2.3 

3.25 

20 

4.6 

3.60 

30 

6.9 

3.75 

40 

9.2 

3.70 

50 

11.5 

3.65 


°C 


(cm^/sec) 




For 30% KQH solution, the diffusivity of the electrolyte is 3.75 x 10 

2 

cm /sec. With the available data on diffusivity and the estimated 
values of the hydration number, the concentration gradients of "free" 
water at different KOH concentrations can be calculated as a function 
of current density. The results of "free" water concentration gradient 
are summarized in Table 5.5, 


Table 5.5 

v"Free" Water Concentration Gradients 


Hydration 
No . of 

Hydroxyl Ion 

KOH 

Concentration 


s_ 

c=2 . 3 

. A. 6 

' 6.9 

9.2 

11.5 

1 

0. 64i 

0. 58i 

0. 55i 

0.561 

0. 57i 

2 

0 . 96i 

0. 86i 

0. 83i 

0.8A1 

- 

3 

1 . 27i 

1 . 15i 

1.101 

- 

- 

A 

1 . 59i 

l.AAi 

- 

- 

- 

5 

1 . 91i 

1. 73i 

— 

- 

- 


For the concentration gradient of the KOH solution, we have from 


Equation 5.2-21 


N 



dine 
o 

dine 



it 

IT 


From Faraday ' s Law , 
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Hence the concentration gradient of KOH solution in the electrolyte 
matrix is 


dc 

dx 


(1 - t ) i 


v D 


1 - 


dine 
c 

dine 


(5.5-2) 


The term 


1 - 


dine 
c 

dine 

( 


1 - 


can be shown to be equal to 


dine 

o 


dine 


(5.5-3) 


c V 
o o 


where V is the partial molar volume of the solvent and is given by 
o 


M 


V = 
o 


p - c 


d£ 

dc 

3 


(5.5-4) 


The density of KOH is 1.1884 gm/cm for 20% KOH solution and 1.2905 
gm/cm 3 for 30% KOH solution. If we assume that the density is a linear 
function of the concentration, we get 


V 

o 


1.2905 gm/cc 


18 gm/mole 


6.9 moles/1 


1.2905 - 1.1884 


2.3 


) 


gin/ cc 
mole/1 


= 18.23 ce/mole 


Hence 


( 


dine 
o 

dine 


1 

50.2 mole/1 x 


18.23 gm/cc 
1000 cell 


1.093 


Thus, the concentration gradient of the KOH solution in the electrolyte 


matrix is given by 
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dc _ 0. 089i 

3.75 x 10 ^ cm^/sec x 1.093 x 96500 amp-sec/mole 

= 0.023i mole/l/cm. 

The concentration gradients of KOH solution at different concentrations 
can be estimated by the same procedure and the results are summarized 
in Table 5.6 


Table 5.6 


KOH 

Concentration 

Gradients in the 

Electrolyte 

Matrix 




D x 10 5 

dine 

- o 

dc 


KOH 

t + 

(cm^/sec) 

dine 

dx 

10 


0.142 

3.25 

0.868 

0.052i 

20 


0.112 

3.60 

0.965 

0 . 033i 

30 


0.089 

3.75 

1.093 

0.023i 

40 


0.070 

3.70 

1.387 

0.0141 


From Table 5.6 we can see that, for the estimated values of the 

transference and hydration numbers, a current density of one ampere 
2 

per cm , and a mean KOH concentration of 30 Wt . % (6.9 N) , the KOH 

concentration gradient in the electrolyte matrix is 0.023 g mole/l/cm. 

Since the electrolyte matrix is normally less than one millimeter in 

thickness, it is clear that the KOH concentration difference between 

the anode and the cathode even at the high current density of one 
2 

ampere per cm appears to be almost negligible. Moreover, it is clear 
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that diffusion and electromigration are entirely adequate to support 
very high current densities. 

5. 6 Discussion of Simplifying Assumptions 

In calculating the concentration gradients, we have assumed that 
convection is absent. From the calculations obtained, the concentration 
gradient of the electrolyte across the matrix is very small. The 
presence of convection would only reduce the concentration gradient 
further, therefore, the assumption of no convection has not led to 
significant errors in estimating the free water concentration gradient. 
In general, a temperature gradient in the electrolyte matrix could 
lead to a variation in composition which would promote convection. 

The temperature gradient can also contribute to changes in transport 
properties and lead to thermal diffusion. Since the transport rates 
are large and convection has been shown to be of negligible importance, 
and since thermal diffusion effects are usually of only secondary 
importance, the assumption of constant temperature seems reasonable. 



6 . 


Effectiveness Factors In Fuel Cell Components 


Porous diffusion electrodes are of extremely complex structure, 
and the sizes of the pores in these electrodes range from 2 ym 
(2 microns) down to as low as about 2 nm (20 A). Electrolyte matrices 
are less complex, but readily accessible to exact analysis. In such 
cases, the rate of solute diffusion within these pores may become 
appreciably less than expected when the solute molecular size becomes 
significant with respect to the pore size. The ratio of the effective 
diffusivity to the bulk diffusivity of any solute in a porous medium is 
termed the effectiveness factor, £, This factor is a function of the 
porosity and tortuosity of the porous medium, of surface effects, 
and of possible hindrances to diffusion owing to the size of the 
diffusing molecules approaching that of the pores. The porosity and 
tortuosity are macroscopic effects; they correct for the reduction in 
the area available for diffusion, and the increase in length of the 
diffusion path, respectively. The other effects are highly complex 
functions of the molecular properties of the solute and porous medium, 
and these are not yet well understood. Beck and Schultz (30) as well as 
Satterfield et al . (31), found that the diffusivity of solutes was 
reduced when the solute diameter approached the pore diameter in their 
studies of diffusion through well-defined, homogeneous porous media. On 
the other hand, Anderson and Quinn (32), and also Wade (34), found no 
alteration of the transport properties in their studies of the influence 
of structure of the porous medium on transport processes in pores with 
sizes comparable to the liquid molecules. 
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It is evident from these contradictory results, and the lack of 
understanding of the influence of other parameters on the effectiveness 
factor in porous electrodes, that experimental measurements of effective- 
ness factor are necessary for a fuller understanding of transport 
processes in electrodes, and it is possible that structural influences 
may be important in the transport processes occurring in electrolyte 
matrices. The sections following describe a proposed method of 
evaluating effectiveness factors experimentally and a solution to the 
mathematical model is given. 

6 . 1 Development of the Model 

The proposed experiment involves measuring both the bulk and 
hindered dif fusivities using a modified stagnant microelectrode . With 
this method only a single measurement is needed for measuring the ratio 
of bulk and hindered dif fusivities . 

A schematic diagram of the diffusion cell is shown below. 



I is a region where bulk diffusion prevails, II is the sample whose 
effectiveness factor is to be measured, and III is another region of 
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bulk diffusion. For such a system, the flux at the electrode can be 
obtained by solving the following set of differential equations, with 
the accompanying initial and boundary conditions, for t 0 


rv 2 , 

3 c x 
3x 2 

1 r— * 
i 

3c| 

3t 

0 

-(^ + l 2 ) < X < -l 2 

(6.1-1) 

3 c 2 
3x 2 

1 

D 2 

3c 2 

3t 

0 

-l 2 < x < 0 

(6.1-2) 

a 2 « 
3 C 3 

3x 2 

L- 

D 3 

3c 3 

3t 

0 

0 < x < l 3 

(6.1-3) 


with the following initial and boundary conditions: 

Initial condition : t <_ 0, c^ = c^ = c^ = c^q for all x. 


Boundary conditions : t > 0 

cj - 0 at x = + 1^) 


c f — c 1 
1 2 


3c' 3c' 

°1 37 " = D 1 "37 


3c' 3c' 

c 2 - C 3> D 2 "77 " D 3 “37 * 


at x = -t r 


at x - 0 


c 1 = c 
3 10 


at x = t 


where c^, c^, c^ are the concentrations in regions I, II, and III, 
respectively. (They are actual concentrations). , and are 

corrected diffusion coefficients in regions I, II, and III, respectively, 
i.e., correction is made for the fact that the flux is based on total 


geometrical area. 



65 


The set of equations above can be solved by Laplace Transformation 
(complete solution is given in Appendix 4), after solving for c^, and 
given that 


3c, 

I_ . = n FAD - 5 -^ 

lim 1 dx 


x - - (Z^ + Z^) 


(6.1-4) 


We obtain the following expression for f^im 1 


+ l 


llm fl . ^3 m =0 

D 1 D 2 D 3 


. 2 

1 t 6 
m m 


m 


(6.1-5) 


where 0 = K_ sin A £ (cos K-A f. 0 tanK 0 A £ + sin K..A ■£„) 

m 1 ml 1 m 2 2m3 1 m 2 

- cos A Z- (K^sin K. A £ tan K_A £_ - K-K- cos K HJ 
mil lm2 2m3 12 , lm2 

P = K 1 (.f 1 + £ ) tan K 0 A £ (sin A Z, cos K A L 
m 11 2 2 m 3 ml lm2 

+ K- cos A Z-. sin K.A Z 0 ) + K 0 (£- + £„) sin A L sin O L 

1 ml • 1 m 2 2 1 2 ml 1 m 2 


( K 2^2 

K„ 


12 1 


+ K„K cos A £, cos K,A Z 


WUU XV- /V v-rt 

ml 1 m 2 


V_3 

K„ 


( K- sin A Z . sin K A - cos X cos K.A Z n \ 
I _L ml 1 m 2 m 1 1 m 2 

2 

cos K_X £„ 

2 m 3 





K 1 D. 


D, 


k 9 =Vtt 
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is given by the equation 



tan K-A £ 0 + K_ tan A K Z. + K-K tan A 

Xm2 1 m 2 3 12 ml 

tan A «£.. tan K-A tan K_A £_ = 0 
ml lm2 2m3 


6. 2 Proposed Equipment and Procedure 


As noted earlier, the proposed method of measuring effectiveness 
factors is based on measuring the limiting current with a system 
approximating a stagnant microelectrode. However, the sample whose 
effectiveness factor is being evaluated will be sandwiched between two 
porous disks in which bulk diffusion prevails, i.e., the pores are so 
large in these regions that there is no hindrance to diffusion in the 
pore. 

In practice we propose, as is often done, to measure the limiting 
current for an electroactive species whose bulk diffusivity is already 
known accurately, and for this purpose we have chosen potassium ferro- 
cyanide. We shall make the measurement using 0.05 M K^Fe(CN)^ in 
0.1 N KC1. 


A somewhat enlarged cross-sectional diagram of the cell is shown 
below. 


Nylon Electrode Body 




Porous discs I and III and the porous sample whose effectiveness 
factor is being measured are fitted into the inside of the cell, 
and" by means of suction the potassium f errocyanide solution is drawn 
into the cell. A platinum electrode is then pressed on the surface of 
disc I by means of a screw-on cap. The experiment can then be started 
by applying a voltage across the diffusion electrode and a counter 
electrode while both are immersed in the potassium ferrocyanide 
solution. A plot of current vs time is generated in the polarograph, 
and the effectiveness factor of the sample is then calculated from the 
experimental curve by fitting Equation (6.1-5), using a trial and error 
procedure on the computer to select appropriate values for and noting 
that (recalling that is already known) 
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APPENDIX 1 


Derivation of Flux Equations (Equations 5.2-11 and 12) 


From Faraday's Law, we have 


i * F l z.N. 
i 1 1 


= F y z .c.v. (A. 1-1) 

“ ixi 
i 

Since the solution, as a whole, is electrically neutral, we have 


y z . c . - o 

H ii 
i 

Thus Equation A. 1-1 can be written as 

i *= F y z.c.v. - F 7 (z.c.)v 
“111 .110 

i l 

= F y z.c. (v . - v ) 

“11 1 o 

i 

Therefore, for a 1-1 electrolyte such as KOH, we have 

i = Fz, c, (v, - v) + Fzc (v - v) 

+ ++ o — — — o 


(A. 1-2) 


(A. 1-3) 


(A. 1-4) 


From Equations 5.2-6 and 7, we have 


c + Vu + = RT 


c - c + 


(V - V ) + -^r 

C 1 0+ C T +- 


(v . - v o ) 


(5.2-6) 


c c c c 

cVu » rt | -r~- (V. - V ) + 


C A- ° 

Adding the last two equations, we get 


Vfc <v » ' 


(5.2-7) 
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c , Vy , + c Vy = 

I i *” ” 


RTc, c RTc c 

““7T (V ~ V ) + — ^ 

C T^+ 0 + C T^- 


(v o - 


v ) 


(A. 1-5) 


The electrochemical potential of the electrolyte, ]i £ j is given by 


u = V ii, + V 11 
e t -r — — * 


(A. 1-6) 


The stoichiometric concentration of the electrolyte can be written as 


c *» 


It = 1 

V. V 


(A. 1-7) 


Hence 


V+ + c -^- 


= c( v + H + + vjl_) 


cVy 


Therefore, we get 


RTc.c RTc c 

c Vp = - (v - v ) + (v - v ) 

e Cn&r, o + c m ' o 


'T o+ 


T o- 


(A.l-8) 


Eliminating v + from Equations A. 1-4 and 8, we get 


RTc, c i 

+ Fz c.cVy 

C T^ + + ' 


RTc c RTc c 

— /v-° Fz c ( v - v ) + — 7 p - ° Fz c (v - v ) 
C T o+ “ - " ° C T ^1 + + o 


RT Co c + c_F i z^ 


z \ 


8" j 

o- J 


(v o - V_) 


X 
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RTc c.c F 

O + — 


z, 3-, - z fr* 

+ 0+ ~ Q“ 


( v o - vj 


Therefore 


iz ccV]i z. 

/ * _ - o- T e 4- o+ o- 

C “ V o V “ z F(z /^ - z ^ ) RTc (z, - z (^ ) 

-+o+ - o- o+o+ - o- 


Substituting Equations 5.2-13 and 14 into Equations A. 1-10 we get 

it_ C T C ^o Z + n 

c _(v o _ v _) = _ _ + ____ - A" 

O + - 

Since the solution is electrically neutral 


z. c, + z c = 0 
+ + — — 


z + c. 


E + - E. 


c. + c 
+ 


V , C + V c 

+ — 


v t + v 

+ 


V 


(A. 1-12) 
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Substituting Equation A. 1-12 into Equation A. 1-11, we get 


c (v - v ) = 
T- o - 


it c cVu V 

- C T p e - ^ 

z F + VRTc 


(A. 1-13) 


Hence 


v_/6e T it_ 

c - v - - - TMr c ^e + 7T + c - v 0 

o 


(5.2-12) 


Similarly, it can be shown that 


V + ^ T 


it 




C + V + “ VRTc ’ z , 

o + 


cVy + — 7 T + c , v 
e z . F +o 


(5.2-11) 



APPENDIX 2 


Derivation of Equation 5.2-16 


The electrochemical potential of the electrolyte is given by 

y = vRTln (cf , a^) (A. 2-1) 

e ± ± 

where f is the mean molar activity coefficient of the electrolyte 
0 

and a + is a proportionality constant. The gradient of the electro- 
chemical potential of the electrolyte is then given by 

Vp = VRT f- ln(cf a®) 
e dx ± ± 


Therefore 


yRT 

c f + a® 


f 0 9c 6 3f ± ] 

f ± a ± 37 + ca ±.37- j 


VRT 

cf 


T r 

+ , 


± + ° 3c J 


Vc 


Thus 


cVy 

e 

VRT 


f 


1 + 


c9f 


Vc 



9lnf + 
9 Inc 


) Vc 


^c T 

vRTc cVU e 22 



f 1 + 




dlnf + 
dine 



(A. 2-2) 


(A. 2-3) 


The measured diffusion coefficient of the electrolyte, D, is related 
to fr'by 
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dlny + 

D = dll + VM m) (1 + — ) (A. 2-4) 

o dlnm 

V = v + + V 


Therefore 


m = 


p- I 
j-0 


c .M . 
3 3 


c 


c M 
o o 



f r 1 


r dlny + ' 

d =tr 

; i + (v + + v_) M o c M 
1 0 0 ; 


dlnm 


-Jp 


= jy ~ 

c 

o 


C + C_|_ + c 


1 + 


1 

J 

dlny^ 

dlnm 


1 + 


dlny ± ' 
dlnm 


Hence we get 


vRTc cVP e = 


D 


1 + 


dlnf + ^ 
dine 


1 + 


dlny d 

dlnm 


Vc 


(A. 2-5) 


(A. 2-6) 


(A. 2-7) 


The measured diffusion coefficient of the electrolyte, D, is also 
related to /£^by the following relation 


D “ 


& 


1 + VM m 
o__ 

1 4* Mm 


r dlnf + 

1 + — 

dine 


dlnp 

dine 


(A. 2-8) 


From Equations A. 2-4 and 8 


we get 
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1 + 


dlnf, 

dine 


1 + 


dlny + 

dlnm 


(1 + Hm) (1 - ^ ) 


1 + 


Me 

c M 
o o 


1 - 


dlnp 

dine 


= 1 + 


Me 

c M 
o o 


f 1 + Me 1 i-M 

^ p I dine 
o 


= 1 + 


Me 


c M 
o o 


pdlnp 

p dine 
o 


= 1 - 


pdlnp 

p dine 
o 


Me 

P_ 


p = c M + Me 
o o 


M 1 -I. M " (C M + MC) 

dc c M + Me oc o o 
o o 


3 c 


M 


p o 3c 


o + M 


Multiplying each side of the above equation by c, we get 


3 c 

dlnp = c M o + 

dine p o 3c 


Me 

P 


cp 3 c x , 

o o Me 

c^pSc p 


P hC y, 

o o Me 

pine p 


(A. 2-9) 


(A. 2-10) 


Substituting Equation A. 2-10 into Equation A. 2-9, we get 
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dlnf ± 

^ + dine 

dlny + 

1 

dlnm 


f 

p 1 

t p dine 

\ ° ...... 9. + 

Me 1 

1 

Me 

P o 1 

[ p dine 

P J 

’ P o 


r dine ,, .. ^ 

o Me Me 

dine p p 

r o o 


dine 
o 

dine 


Substituting Equation A. 2-11 into Equation A. 2-7, we thus get 




VRTc 





dine 'i 

o 

dine ) 


Vc 


(A. 2-11) 


(5.2-16) 



APPENDIX 3 


Sample Calculation of "Free" Water Concentration 


The amount of "free” water for a given electrolytic solution is 
equal to the portion of water that is not tightly bound in the primary 
hydration sheaths of the ion. Hence the concentration of "free” water 
is the concentration difference before and after hydration of ions. 

For 30% KOH solution, the concentration of KOH and water before 
hydration are as follows: 


Concentration of KOH, 


c 


30 gm 

100 gm soln 


56.1 


gm 

mole 


= 6.9 gm/1 


The molar concentration of water 

- 30) g m X 1.29 g™ SOl - n ' X 1000 ^ 

100 gm soln cc 1 

18 - g ~- 
mole 

= 50.2 mole/1 

If the hydration number of the potassium ion and that of the hydroxyl 
ion are assumed to be equal, then the concentration of "free 11 water is 
given by 

C fw °tw °bw 


= c 


tw 


2sc 
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where c is the concentration of total water, that is, the concentration 
of water before hydration; c^ w is the concentration of bound water. 

If the hydration number of the ions is equal to one, then 

C fw = 50,2 mole / i ~ (2)(1)(6.9) mole/1 

= 36.4 tnole/l 


Values of the concentration of "free" water are summarized in Table 5.2. 



APPENDIX 4 


Solution of the Diffusion Equation 
For the Effectiveness Factor Cell 

The set of equations to be solved is, for t > 0 

0 ~(£ 1 + < x < 

0 -l < x < 0 

0 0 < x < l 3 

where the concentrations are based on the actual concentration in the 
pores, and the diffusion coefficients are those corrected for the 
differences in diffusion surface area. 

The boundary and initial conditions for the above system are; 

Initial Conditions: t 0 

c » = r * — r * = r 

C 1 2 3 10 

+ l 2 ) 

at x = - ^ 


Boundary Conditions: 


C 1 = 0 


at x = 


c f = c f 
1 2 


3c* 

*i a^T 


= D 


^2 

2 3x 


3 2 cj 

1 

3c i 

3x 2 

' D 1 

3t 

3 2 c; 

1 


3x ? 

D 2 

3t 

3 c 3 

1 

3c 3 

3x 2 

‘ D 3 

3t 
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c f = c f 
2 3 


at x = ■£. 


D, 


3c 2 


= D, 


3c 3 


2 3x 3 3x 

First, we make the transformation 


c i ~ c i ~ c io 


C„ — — C- ~ 

2 2 10 


c = c' - c 

3 3 10 


The above equations are transformed into 


r.2 

3 c x 
3x 2 

1 

" D i 

3 

3t 

0 

+ l 2 ) < x < -l 2 

(A. 4-1) 

~2 

3 C 2 
3x 2 

i 

D 2 

3c 2 

3t 

0 

-l 2 < x < 0 

. (A. 4-2) 

3 C 3 
3x 2 

1 

°3 

3c 3 

3t 

0 

0 < x < 

(A. 4-3) 


Initial Condition 


t < 0 


C 1 “ c 2 = C 3 = 0 


t > 0 


C 1 = C ls (= ~ c 10 ) 


at x 


= -U 1 + l 2 ) (A. 4-4) 


Boundary Condition: 



81 


" V 

D 1 

9c^ 

3x 

D 2 

3c 2 

3x 

x = -t 2 

= C 3> 

D 2 

3c 2 

3x 

’ D 3 

Sc 3 

3x 

X = 0 


To solve Equations A. 4-1, 2, and 3, with the given Initial 
Conditions and Boundary Conditions, they are Laplace Transformed, 
and solved to give 


c i = 


'Is 


cosh + x) + A sinh q^ + x ) (A. 4^8) 


= B cosh q^ (-^2 + x) + C sinh q 0 (£ Q + x) 


2 K 2 


(A. 4-9) 


= E (sinh q^x - tanh q^-t^ cosh q^x) 


(A. 4-10) 


where 

«1 - ^ 

q 2 = /p/D 2 


q 3 = 

p = the parameter for the Laplace Transformation. 


A, B, C and D are constants to be evaluated from the Boundary 
Conditions (A. 4-5 and 6). It should be noted that conditions A, 4-4 and 7 
have been used to obtain Equations A. 4-8 through 10. 

Substitution of Equations A. 4-8 through 10 into the Boundary 


Conditions give 
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c 

— — = cosh q + A sinh q £ = B 

p 1 -L x 

(A. 4-11) 

C 1 

» D 1 q 1 sinh q^ + Aq^D 2 cosh q^ = cq 2 D 2 

(A. 4-12) 

B cosh q^ + c sinh q^ = - E tanh q^ 

(A. 4-13) 

B q 2 D 2 sinh q 2 -£- 2 ^ q 2^2 q 2^2 — ^^3 q 3 

(A. 4-14) 


Solving Equations A. 4—11 through 14 for A gives 

A = fiE { m/N } (A. 4-15) 

P 

where M = sinh (^ 3^3 sinh <^2 + D 2 q 2~- COSh q 2'^2' ) 

+ cosh q^£^ ^ D 3 q 3 cos ^ ^ 2^2 + D 2 q 2 sin ^ q 2^2 tanh q 3^V 

N = D 1 q 1 cosh q^ (D^ sinh q 2 ^2 + D 2 q 2 C0Sh q 2'^2 tanh q 3 ^' 3 ‘ ) 

+ D q sinh q^ (D^ cosh + D 2 q 2 sinh q^ tanh q^) 


Only A was solved for because we are interest ed ' in c^ which contains A 
as the sole constant to be evaluated. 

Substituting Equation A. 4-13 to Equation A. 4-8 yields 



CA.4-16) 


where K = (D^ sinh q^ + cosh q^ tanh q^) cosh q^l^lfx) 

- Vq (D 3 q 3 cosh q^ + D 2 q 2 sinh q^ tanh q^) sinh q-^-^+x) 
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L = ^ D 3 q 3 sinh ^2^2 + D 2 q 2 cosh q 2^'2 tanh q 3'^3^ costl q l^l 

+ D 2 q 2 (D 3 q 3 cosh q^ + sinh q^ tanh q^) sinh q^ 

can be obtained by means of the Inversion Theorem. We see that 
Equation A. 4-16 gives a simple pole, and zeros of "L" give simple 
poles of the integrand of Equation A.4VL6 in the inversion equation. 


— - Res (0) + l Res (-D n X 2 ) 

, L . 1 m 

Is m=l 


We see that Res (0) = 0/0; therefore, we have to apply L’Hopital’s 
rule, which, after simplification, gives 

D 1 D 3 £ 2 + D 1 D 2^3 ~ D 2 D 3 (£ 2 + x) - 
Res (0) " Dl D 3 £ 2 + Dl D 2 £ 3 + D 2 D 3 £ x 


r 

Res (-D.X ) = 3 


(By using L'Hopital's rule) 


OL/3p) _ D 


1 m 


Letting = ^2^1 * ^2 ^3^1 * 


p = -D 1 X 
1 m 


2 

-iX t 

K(-D^X 2 ) = ie m ( /D 1 D 2 cos 1(2^ + l 2 + x) cos K^Xt tan K 2 Xt 3 

D 2 D 3 

+ v'dTdT cos X(2 + l n + x) sin K,X^ n - sin X(£ 0 + x) cos K n X£ 0 

12 12 -12 2 1 2 


+ D 2 sin XK.^ sin X(i- 2 + x) tan K 2 ^3^ 
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dL 

dp 


+ 


+ 


where 


p=-D 


_ 1 
2 “ 2i 


1 m 


sin AL^ cosAK^^ tan ^ K 2^3 
“^2 cosAL^ sin tan ^^2^3 


cos AL^ cos 


K 2 3 


cos K^AZ^ 


- + ^ 2 ^ Sin ^1 Sin ^ K 1^2 


K 


~ Jt^ cos A cos AK^^ + C0S ^1 C ° S ^1^2 


K 2 L^ cos Xt^ sin AK^-^ tan ^ K 2^3 


K ^-^2 s: *- n COS ^ 1^2 tan ^^ 2^3 


a 

K 2 3 


sin AL^ sin 
cos 2 AK 2 £ 3 


'Is 


D..D l + D D L - D D (L + x) 

w, + w 3 - Wi - 


I 2 e 


D-A t 9 
1 m m 


m 


m 


0 = {K. cos A (2 1. + £ 0 + x)(cos K A l tan K A L 

m 1 m l 2 lm2 zmj 


+ sin K- A Z n ) + sin A (L 0 + x)(K 1 sin K.A L tan KJ L 
1 m 2 m2 1 lm2 zmo 


- K.iC cos KA £-)} 
12 m 2 


P = D,A J K- (L, + L 0 ) sin A L- cos A K L tan A K L 
m 1 m < 1 1 2 ml m 1 2 m23 


+ K 2 (L- 4- £ 0 ) cos A £- sin A K L tan AK l 
l v 1 2 ml m 1 2 23 

( K 2^2 

+ + V sin sin X^ - 


K. + K lVl 


(K.. sm A L n sin A - cos A L- cos A K L) 
, „ , „ , 13 M ml m ml ml2 

cos A L- cos A K £. 0 + ~ „ — o 

ml m 1 2 K„ 


cos A K,JL 
m 2 3 
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where A^ are given by solution of the equation 


K 2 tan K x A£ 2 + K ± tan K 2 A£ 3 + i^K tanA^ 


+ tan K^A£ 2 tan ^A^ 3 tanA£^ = 0 


since 


hi™ - nFAD i w 


X = ~ l l - l 2 


Finally, we have 


'lim 


c 10 nFA 


t l ? l. 

\ + \ + S 


iA t e 

. o V m 111 
+ 2 i e p- 

m m 



LIST OF REFERENCES 


1. ^Feldkamp, K. , Chemie. Ing. Tech,, 41 , 1181 (1969). 

2. International Critical Tables, 1st Ed., IV, 466 (1928). 

3. Faust, Von 0., Zeit. Anorg. Chem., 160 , 373 (1927). 

4. Walker, R. D. , Fifteenth Semi-Annual Report, September, 1973, 

NASA Research Grant NGR 10-005-022. 

5. Kordesch, K, V. , From Electrocatalysis to Fuel Cells . Washington: 
Washington University Press, (1972), p. 157. 

6. Miller, M. L., and Fornasar, H. J. , J. Electrochem. Soc., 115 
330 (1968). 

7. Lundquist, J. T. , and Vogel, W. M. , Ibid., 116 , 1066 (1969). 

8. Bockris, J. 0’M, Quart. Rev. (London) _3, 173 (1949). 

9. Hirschf elder , J. 0., Curtiss, C. F. , and Bird, R. B. , Molecular 

Theory of Gases and Liquids , New York: John Wiley & Sons, Inc., 

(1954) p. 718. 

10. Bird, R. B. , Stewart, W. E. , and Lightfoot, E. N., Transport 

Phenomena . New York: John Wiley & Sons, Inc., (1960) p. 570. 

11. Maclnnes, D. A., and Dole, M. , J. Am, Chem. Soc., 53^, 1357 (1931). 

12. Jones, G., and Bradshaw, B, C. , Ibid., S4, 138 (1932). 

13. Maclnnes, D. A., and Longsworth, L. G. , Chem. Rev., 11 , 171 (1932) 

14. Jones, G, , and Dole, M. , J. Am. Chem. Soc., JLL, 1073 (1929). 

15. Longsworth, L. G. , Ibid., 54 , 2741 (1932). 

16. Owen, B. B,, Ibid., JT7, 2441 (1935). 

17. Merenkov, P. T. , Uzbeksh Khim. Zn. 50, (1961). 

18. Knobel, N., et al . , J. Am. Chem. Soc., 4_5, 77 (1923). 

19. Ulich, H., Z. Elektrochem. , 36, 497 (1930). 

20. Darmois, E. , J. Phys. Radium, _8, 117 (1942). 


86 



87 


21. Passynski, A. , Acta Physiocochim. , _8, 385 (1938) . 

22. Bockris, J. O'M. , and Saluja, P. P. S., J. Phys . Chem., 76 , 

2140 (1972) . 

23. Hindman, J. C. , J. Chem. Phys., 36_, 1000 (1962). 

24. Creekmore, R. W. , and Reilly, C. N. , J. Phys. Chem., 7^3, 1563 (1969). 

25. Brady, G. W. , and Krause, J. T. , J*. Chem. Phys., 21 _ , 304 (1962). 

26. Swinehart, J. H., and Taube, H. , J. Phys. Chem., 31_, 1579 (1962). 

27. Jackson, J. A., and Lemons, J. F. , and Taube, II., J. Chem. Phys., 

32, 553 (1960). 

28. Isemura, Toshizo, and Goto, Sachlo, Bulletin of the Chemical 
Society of Japan, _3_7, 1690 (1964) . 

29. Bhatia, R. N. , M. S. E. Thesis, University of Florida (1960). 

30. Beck, R. E. , and J. S. Schultz, Science, 170 , 1302 (1970). 

31. Satterfield, C. N., C. K, Colton and W. H. Pitcher, A. I. Ch. E. J, 
19, 628 (1973). 

32. Anderson, J. L. and J. A. Quinn, Trans. Farad. Soc. I. 68 , 

744 (1972). 

33. Wade, W. H. , Soc. Pet. Eng. of A. I. M. E. SPE 4364 (1973). 

34. Walker, R. D. , Sixth Semi-Annual Report, February, 1969, NASA 
Research Grant MGR 10-005-022. 


x 



